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Abstract. An alternative method for solving the fractional kinetic equations 
solved earlier by Haubold and Mathai (2000) and Saxena et al. (2002, 2004a, 
2004b) is recently given by Saxena and Kalla (2007). This method can also 
be applied in solving more general fractional kinetic equations than the ones 
solved by the aforesaid authors. In view of the usefulness and importance of 
the kinetic equation in certain physical problems governing reaction-diffusion 
in complex systems and anomalous diffusion, the authors present an alternative 
simple method for deriving the solution of the generalized forms of the fractional 
kinetic equations solved by the aforesaid authors and Nonnenmacher and Met- 
zler (1995). The method depends on the use of the Riemann-Liouville fractional 
calculus operators. It has been shown by the application of Riemann-Liouville 
fractional integral operator and its interesting properties, that the solution of the 
given fractional kinetic equation can be obtained in a straight-forward manner. 
This method does not make use of the Laplace transform. 

1 Introduction 

The paper deals with the essential problem related to applications of Mittag- 
Leffler function and Riemann-Liouville fractional calculus operators to fractional 
order kinetic equations arising in modeling physical phenomena, governing diffu- 
sion in porous media and relaxation processes. As such it reveals the important 
role of these tools in applications of fractional calculus. The results are inter- 
esting and useful for wide range of applied scientists dealing with fractional 
order differential and fractional order integral equations. In a series of papers 
the authors have demonstrated the use of integral transforms in the solution 
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of certain fractional kinetic equations (2002, 2004a 2004b), reaction-diffusion 
equations (2006a, 2006b), and fractional differential equations governing non- 
linear waves (2006c, 2006d ). In the present paper it is shown by the application 
of Riemann-Liouville fractional calculus operators and its interesting properties 
that the given fractional kinetic equations can be easily solved. Fractional ki- 
netic equations are studied by Zaslavsky (1994), Saichev and Zaslavsky (1997), 
Glocckle and Nonnenmachcr (1991), and Saxena, Mathai and Haubold (2002, 
2004a, 2004b) due to their importance in the solution of certain applied problems 
governing reaction and relaxation in complex systems and anomalous diffusion. 
The use of fractional kinetic equations in many problems arising in science and 
engineering can be found in the monographs by Podlubny (1999), Hilfer (2000), 
and Kilbas, Srivastava and Trujillo (2006) and the various papers given therein. 
The Mittag-Leffler functions naturally occur as a solution of fractional order dif- 
ferential equation or a fractional order integral equation. Mittag-Leffler (1903) 
defined this function, known as Mittag-Leffler function in the literature, in terms 
of the power series 

E -W := E r( afc + l) ; (« e ^ e ( a )>°)- 
This function is generalized by Wiman (1905) in the form 

T(ak + p) 



E a ,p{z) := J2 WTTWY (a ' £ G °> Re{a) > °' Re ^ > 0) ' (2) 



According to Dzherbashyan (1966, p. 118), both the functions defined by the 
equations (1) and (2) are entire functions of order p = l/a and type a = 
1. A comprehensive detailed account of these functions is available from the 
monographs of Erdelyi, Magnus, Oberhettinger and Tricomi (1995, Chapter 
18) and Dzherbashyan (1966, Chapter 2). The Riemann-Liouville operators of 
fractional calculus are defined in the books by Miller and Ross (1993), Oldham 
and Spanicr (1974), Podlubny (1999) and Kilbas, Srivastava and Trujillo (2006) 
as 

1 /•* 

a Dt V N(t) := — - (t- u)"- 1 N(u)du, Re{v) > 0,t > a (3) 

^ V^J J a 

with a D° t N(t) = N(t), and 

a D?N(t) := —( aD r n N(t)),Re(n) >0,n- f i>0. (4) 
By virtue of the definitions (3), it is not difficult to show that 

a D^(t a)^ 1 = p^y(i - ay + -\ (5) 
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where Re(u) > 0, Re{p) > 0; t > a. Also from (Podlubny, 1999, p.72, eq.(2.117)), 
we have 

a D»(t a)"- 1 = - a)"-"" 1 , (6) 

where Re{v) > 0,Re(p) > 0, t > a. When p = 1 (6) reduces to an interesting 
formula 

a Dt v l = \ - (t - a)~\t > a; v ± 1, 2 . . . (7) 

which is a remarkable result in the theory of fractional calculus and indicates 
that the fractional derivative of a constant is not zero. 

We now proceed to derive and solve the fractional kinetic equations in the 
next section. 



2 Derivation of the fractional kinetic equation 
and its solution 

If we integrate the standard kinetic equation 

j t N i (t) = -c i N i (t),(c i >0) (8) 

we obtain (Haubold and Mathai, 2000, p. 58) 

N(t) - N(a) = -a aD^NM, (9) 

where aD^ 1 is the standard Riemann integral operator. Here we recall that in 
the original paper of Haubold and Mathai (2000), the number density of species, 
Ni — Ni(t) is a function of time. Further we assume that Ni(t = a) = N a is 
the number density of species i at time t = a. If we drop the index i in (9) and 
generalize it, we arrive at the fractional kinetic equation 

N(t) -N a = -c v a D^N(t) (10) 

Solution of (10). If we multiply both sides of (10) by {-c u ) m a D^ mv \ we obtain 

{-c v ) m a Dt m »N{t)-{-c v )(-c v ) m a Dt mv ~ v N(t) = (-c v ) m a Dl mv N a . (11) 

Now summing up both sides of (2.4) for m from to oo, it yields 



{-c v ) m a D- mv N{t) - {-c v ) m+1 a Dt m "- v N{t) 

771=0 777 — 

OO 

= N a (-cT" aA" m "l, (12) 



m=0 



3 



which on using the formula (5) yields 

oo 

N(t) = JV a ^(-O m [(t-a) m 7r(mi/ + l)] (13) 
= N a E v [-c u {t-a) v ],t> a (14) 
Thus we arrive at the following theorem: 

Theorem 1. If Re{v) > 0,i?e(c) > then there exists the unique solution of 
the integral equation 

N(t)-N a = -<? a Dt v (t), (15) 

given by 

N(t) = N a E v (-c v {t - a)"), t > a (16) 

with the Mittag-LefHer function defined by (1). 

When a — > 0, (16) reduces to the following result given by Haubold and 
Mathai (2000, p.63): 

Corollary 1.1. If, Re(c) > then the unique solution of the integral equation 

N(t) - N = -c%£>r (*), (17) 

is given by 

N(t) = N E u {-c v t v ). (18) 
Remark 1. If we apply the operator a D^ from the left to (10) and make use 
of (7), we obtain the fractional differential equation 

a D?N(t) - N a ^ °) " = -c v N{t), t>a (19) 

whose solution is also given by (16). When a tends to zero in (16), it reduces 
to one obtained by Nonnenmacher and Metzlcr (1995, p. 156) for the fractional 
relaxation equation with c replaced by 1 / c. 

Remark 2. The method adopted in deriving the solution of fractional kinetic 
equation (8) is similar to that used by Al-Saqabi and Tuan (1996) for solving 
differ integral equations. 

3 Theorem 2. 

If min {Re{v) , Re(^i)} > 0, Re(c) > 0, then there exists the unique solution of 
the integral equation 
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given by 

N(t) = N a T(fx)(t - a)"- 1 E u ^{-c v {t - a) v ),t > a (21) 

where E v ^(t)\s the generalized Mittag-LcfHcr function defined by (2). 

Solution of (20). If we multiply both sides of (20) by (-c v ) m a Vt ' mv , we 
obtain 

{-c v ) m a D^ mv N{t) - {-c v ){-c v ) m a Di rnv - u N(t) = N a (-c u ) m aD^t^ 1 . 

(22) 

Now summing up both sides of (22) for m from to oo, it yields 

oo oo oo 

m— m — m—0 

(23) 

which on using the formula (5) gives 

oo 

N(t) = 7V a r( M )^(-O m [(t-«) m 77(^ + M )] (24) 
- N a E v ^[-c v {t-ay\,t>a. (25) 

This completes the proof of Theorem 2. 

For a = 0, (25) reduces to the following result given by Saxena, Mathai and 
Haubold (2002, p.283). 

Corollary 2.1. If min {Re(v), Re(fi)} > 0, R{c) > then the solution of the 
integral equation 

N(t) - A^" 1 - -c" oDt"{t) (26) 

' is given by 

N(t) = NoTifi)^- 1 E^{-c v t»), (27) 
where E u ^(t)is the generalized Mittag-Lcfner function defined by (2). 
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